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1. Introduction 

Superstring theories and M-theory can be consistently realized on various classical solutions 
of supergravity theories. However, it is difficult to concretely study theories on these 
backgrounds because of complicated interaction terms. In the recent progress of string 
theory, plane-wave backgrounds are focused upon. The type IIB pp-wave solution was 
found [Ql and it was shown that this background can be described as the approximation 
of the ^4^5' space geometry around a certain null geodesies [§, ^] via the Penrose limit Q . 
The superstring theory on this background |5|, ^ is exactly solvable in the Green-Schwarz 
formulation, in spite of the presence of the Ramond-Ramond fluxes. Moreover, this string 
theory was utilized for the study of the AdS /CFT correspondence at the stringy level Q 
beyond the supergravity analysis. 

On the other hand, the matrix model on the eleven-dimensional pp-wave background 
was proposed by Berenstein, Maldacena and Nastase 0, which is closely related to a 



supermembrane theory on this background M, H via the matrix regularization |10|. This pp- 



wave matrix model includes mass terms and the Myers term [11|, and hence one can expect 
interesting physics intrinsic to the pp-wave case. Because of mass terms, all of the flat 
directions are completely lifted up. So it might be expected that a single supermembrane 
should be stabilized. In addition, the presence of the Myers term leads to many classical 
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solutions living only on the pp-wave. These solutions are fuzzy sphere type ones. In 
particular, the fuzzy sphere solution with zero energy exists, and it can appear in the 
classical vacuum of the system. As the result, the vacuum of this matrix model is very 
enriched. 

Thus, superstring theories and M-theory on plane-wave backgrounds are quite inter- 
esting and so have been intensively studied. However, most of the studies have considered 
the time- independent backgrounds. A recent interest in supergravity, superstring theories 
and M-theory is to include the time- dependent background in our consideration. In general, 
it is much difficult to treat the time-dependent backgrounds and so one can hardly study 
string theories and M-theory on such backgrounds. Under such a circumstance, the time- 
dependent plane-wave background was considered by Papadopoulos, Russo and Tseytlin 
fl^ , and they showed that the string theory on this background is exactly solvable. The 
solvability of this model comes from supersymmetries which are preserved as a character- 
istic of plane-wave type backgrounds. Thus, we can find that the plane-wave backgrounds 
are very available to study the physics on time-dependent backgrounds. 

The main focus of this paper is to study a matrix model on a homogeneous plane- 
wave background and study various features such as supersymmetries, vacuum structure 



and classical solutions in this model. The homogeneous plane- wave background [13| we will 
consider here leads to the time-dependent background through a coordinate transformation. 
Hence the study of the matrix model on this homogeneous background is equivalent to that 
on the time-dependent background. On the other hand, this background can be obtained 
from an M-theory Godel universe* via the Penrose limit and this is also an anti-Mach type 
solution^ |15|. Therefore it can be expected that the study here should be closely connected 
to the M-theory on the Godel universe. 

Our paper is organized as follows: In Section 2, we will briefly review a homogeneous 
plane-wave background which leads to a time-dependent background. In Section 3, the 
action of the matrix model on this background will be proposed. Section 4 is devoted to 
the study of the vacuum structure of the abelian part of the matrix model. In Section 5, 
examples of classical solutions in this model will be presented. We will discuss graviton 
solutions in our model and find fuzzy sphere and hyperbolic type solutions. Section 6 is 
devoted to a conclusion and discussions. In Appendix A, we will discuss a matrix model 
on the general homogeneous plane-wave. Appendix B is devoted to the detail study of the 
energy spectrum. 

2. Time-dependent Homogeneous Plane-wave 

In this section we will briefly review homogeneous plane- wave backgrounds, which are 



Cahen-Wallach space [16| with rotation terms. This type of backgrounds are related to 



time-dependent homogeneous plane-waves via time-dependent coordinate transformations. 



*The Godel universe was originally proposed in the paper 

t M-theory Godel universe can be constructed by lifting up a maximally supersymmetric Godel solution in 
the minimal supergravity in five dimensions jl^. Also, string theories on Godel universe and homogeneous 
plane-wave background are discussed in [|f8| and UM, respectively. 
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These are solutions of the eleven-dimensional supergravity and the family of them is given 
(in the Brinkmann coordinates) by 



ds"^ = -2dx+dx- + fi'^Aijx'x^idx+f + 2fifijx'dx^dx+ + {dx'f , (2.1) 
3! 



F4 = A iijk dx'^^ , (2.2) 



where Aij and fij are constant symmetric and anti-symmetric matrices respectively, and 
^ijk is a completely anti-symmetric third-rank tensor. The nonzero-components of spin 
connection for this background are 

= AijxUx^ + fijdx^ , Lo'^ = -fijdx^ . (2.3) 

The metric ( |2.1[) explicitly describes the time-dependent background if A^j and fij do not 
commute, as mentioned in |jl^. In particular, the above type of plane- wave backgrounds 
is related to the explicitly time-dependent metric 

ds'^ = -2dx+dx- + (e^^^^oe"'^'"^) z'z^{dx+f + {dz'f . (2.4) 



By the use of the coordinate transformation of transverse directions 



z'^w' = U-''n z^ (2.5) 

V / ik 

the metric (2.4) takes the stationary form 

ds^ = -Idx+dx- + ((^o)ii - fikfkj) w'w^{dx+f + {dw'f - 2w' fikdw'' dx+ , (2.6) 



and we see that the time- dependent plane-wave ( P^ ) can be mapped to the stationary 
plane- wave as in ( |2.1| ). 

In this paper, we concentrate on the special case which admits 20 super symmetries 
|2C], where the matrices Aij and fij are given by 

^22=4(P2_1), ^33 = ^44 = -!, /i2 = -P, (2.7) 

and otherwise zero. The (,ijk is written as 

629 = -2, 649 = 2P, 656 = -2(1 - P')^/^ 678 = -2(1 - P2)^/^ (2.8) 

Notably, the above Aij and fij do not commute, and so this background is equivalent to 
the time-dependent background. We can obtain this background from the M-theory Godel 
universe by the use of the Penrose limit Q. The M-theory Godel universe is constructed 
by lifting up a maximally supersymmetric Godel solution of minimal five-dimensional su- 
pergravity [17| (/9 is an arbitrary parameter): 

dsl(. = -{dt + P{rld(j)l + rld(t)l)f + drl + rld(t)l + drl + r^dcpl , (2.9) 

as the product space with the six-dimensional flat space M^: 

9 

dsli = dslc + dz'^ + '^{dx' f , (2.10) 

i=5 
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where this eleven-dimensional background is supported by the four-form field-strength: 



-Pri</>i56 — Fri4)i78 — Fri<j>i9z — Fr2<t>2^& — -^r2</>2 78 — Fr2(l>29z — —2/3 . (2-11) 

The plane-wave solution (2.1) and (2.2) admits 20 supersymmetries (not maximal 32), 
which are inherited from that of the M-theory Godel background (2.7), and its Killing 



spinor was explicitly presented in |2C]. That is, this background admits extra Killing 
spinors satisfying 

[7,^2 + 6^7,0 + 9e\ + 6^i{9, 0] + 18[9, - IU{A,,^^ + fijfjkl^)]e = , 
e{x+) = (^+^'')eo (2.12) 
where we have defined as 

^^^w'' (2-13) 

and eg is a constant spinor. It is worth noting that the parameter P characterizes the null 
geodesic in taking the Penrose limit. The null geodesic is chosen to have a momentum P 
along the six directions transverse to the five-dimensional Godel metric. If we consider the 
P = case, then the null geodesic is in the Godel space in five dimensions. 

3. Action of Matrix Model on a Time-dependent Plane-wave 

Here we shall discuss the matrix model on the general homogeneous plane- wave background 
(anti-Mach type space) whose metric is ( |2.1[) with (|2.7] ) and the four-form flux (^]^) with 
( |2.8| ) is equipped. By generalizing the result of work ||2l|, we can propose the action of 
matrix model in the light-cone gauge as follows: 



S = I dT ^Tr 



- V(l - P^){X^ f - ^2(^3)2 _ ^2(^4)2 

+2fiPX'^DrX^ - 2nPX^DrX'^ + ^i^j'^ {W - ^7^^) ^p 

+ (1 - p2)l/2x[2x5x6] + (1 - P2)l/2^[2^7^8]| 



(3.1) 



W = ^129 _ p^349 ^ (1 _ p2)l/2^256 ^ (^^ _ p2)l/2^278 ^ (3^2) 

where ip is the 16-components 5*0(9) spinor and 7*'s are the 5*0(9) gamma matrices. The r 
is the world-volume time and the covariant derivative Dr is defined as D^-* = dr* —i[uJ, *] ■ 
The gauge connection u describes the area preserving diffeomorphism in terms of membrane 
theory^. It is an easy task to generalize the above action to the background (p.lj) with 



'"We restrict ourselves to the closed membrane case. It is an interesting issue to study open supermem- 
brane theories on a time-dependent background. 
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( |2.2| ), and present the matrix model on the general homogeneous plane- wave (anti-Mach) 
background. This generalization will be discussed in the appendix 

The above matrix model includes the Myers terms due to the presence of the constant 
four-form flux. Hence the fuzzy sphere solutions should appear as classical solutions as in 
the pp-wave matrix model Q. These solutions will be discussed in detail later. 

Here we will discuss supersymmetries of the matrix model on the time-dependent 
background. The 4 dynamical and 16 kinematical supersymmetries are preserved in this 
background, and so the same number of supersymmetries should be preserved in the matrix 
model on this background. 

The transformation law of 4 dynamical supersymmetries is given by 

6,X' = iij^j'e{T), 5,cj = #^e(r), e(r) = et(^+3^^'')"PiP2eo , (3.3) 



where eo is a constant spinor and we have introduced the projection operators: 

Pl = \i^- 75678) , (3.4) 

-P2 = ^ (I + C71956 - S71234) , (3.5) 
where c = (1 — P^)^/^ and s = P. The following properties: 

Pf = Pl, P| = P2, PlP2 = P2Pl. (3.6) 

are satisfied, and Pi and P2 are independent each other. Hence a quarter of dynamical 
supersymmetries (i.e., 4 dynamical supersymmetries) can survive. 

The transformation law of 16 kinematical supersymmetries is represented by 

<5^X^ = V = 0' ^n'i^ = ri{r), r/(T) = e-t(^-^^'')-r?o , (3.7) 

where ryo is a constant spinor. 

We can calculate the superalgebra in our model by using the standard Dirac bracket. 
In fact, our matrix model contains complicated fluxes and the characteristic parameter P 
is turned on, and hence its superalgebra is expected to contain the generators of rotations 
(angular momenta) and nontrivial terms intrinsic to our system. This interesting topic will 



be studied in detail in another place |22 



4. Vacuum Energy of the Abelian Part 

In this section we will consider the abelian part of the matrix model proposed in Section 3. 
This part can be exactly solved and so one can calculate the vacuum energy of this part. 
We will see below that this energy is negative. 
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First, let us consider the P = case. The abehan part of the action is given by 



dT 



• (4.1) 



where we introduced the notations Wq = 7129 + 7256 + 7278 ^-i^d = Ot-X^ 
The dynamical supersymmetry transformation is given by 

S,X' = iV^ye(r) , e(r) = et^o-PiPaco , 



XS - f^X'jiWo - ^X'Woji 



(4.2) 



where Pi and P2 are written as 



Pl = ^il- 75678) , P2 = ^ (I + 71956) • 

On the other hand, the kinematical supersymmetry transformation is described by 

Sr,X' = 0, dr,iP = riiT), r?(r) = e-t^o- . (4.3) 

Here we shall decompose the 16-components spinor as 

^ = ^(++) + ^(+-) + + ipi—) , (4.4) 

where each of the spinors tp^^ \ tp^ and ip^ ^ have only 4 non-trivial compo- 

nents and satisfy the following two chirality conditions: 



75678V'^^*^ = ±V'^^*^ , 71956V'^'^^ = ±V'^'^^ ■ 



(4.5) 



Notably, the matrix Wq commute with 75578 and 71956 , and hence the chirality is preserved 
under the action of Wq (for example, the chirality of WqiIj^'^'^^ is the same with ^(++)). By 
using this decomposition, we can rewrite the fermionic action as 



Sf 



+|i^(++)Tn^(++) _ ^i^(+-)Tn^(+-) 



(4.6) 



where the matrix 11 is defined as 11 = 7255. 

The vacuum energy of bosons comes from the massive directions: X"^, X^ and X^, and 
the net contribution is expressed as 



£;b = I (2 + 1 + 1) = +2// . 



(4.7) 
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The vacuum energy of fermions are evaluated as (for example, by following the procedure 
in the work [p^ ]) 



= -|(l + l + l + 3) = -3^. 
Thus the total vacuum energy is given by 

-E'tot = — • 



(4.8) 



(4.9) 



In conclusion, the vacuum energy of the abelian part is negative^. 

Next, we will consider the P / case. The action of the abelian part is written as 



11^ 



S=-ldr 



(X^y - 4^^(l - P2)(X2)2 _ ^\X^f - ^2(X4^2 



+2fiP{X'^X^ - X'^X^) + ii/j'^i; + ^iij'^ {W - Pj^"^) 



where the expression of W in the P ^ case is 

W = 7^29 

The transformation law of 4 dynamical supersymmetries is given by 
6,X^ = ii;^^'e{T) , e(T) = et(^+3^'7-)rp^p^,^ ^ 



(4.10) 



(4.11) 



(4.12) 



x'-ti + 2ixP{x'Y -x'r) 



where eo is a constant spinor and the projection operators are written as 

-Pi = ^ (I - 75678) , 

-P2 = ^ (I + C71956 - S71234) • 



(4.13) 

(4.14) 



The transformation law of 16 Kinematical supersymmetries is given in (3.7). 

To begin with, we will consider the bosonic part. In the P ^ case the analysis of X^ 
and X"^ directions are nontrivial. We can find that the effect of P is to shift the origin of 
oscillator and to rescale the continuous spectrum of X^, as we will discuss in Appendix. 
Thus, the resulting zero-point energy for bosons is identical with that in the P = case: 



(4.15) 



^The negative vacuum energy of the abelian part of the matrix model on less supersymmetric pp- waves 
is also shown in 1231. 
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Next let us consider the fermionic part and focus upon the fermion mass term: 



i^itp {W - Pj^^) ^ . (4.16) 

We decompose the spinor ^p according to the chirality in terms of the matrices j^^"^^ and 
1956 1234 *i P = case. Remarkably, these matrices commute with W — Pj^"^ . 

We can easily see this property if we rewrite the W — P^^"^ as 

W - Pt'' = 7''' (I + (I - 7'''')(c7'''' + ^7'''')) , 

where we have used the gamma matrices 7*'s satisfy ^^23456789 _ j_ Thus, the chiralities 
of spinor are preserved under the action of — P^^"^ even in the P ^ case. Now we can 
rewrite the fermion mass term as 

where H' = 7-^^^ and the spinors ip~^~^ , ip~^ , ip ^ and 'ip satisfy the following chirality 
conditions with respect to the matrices 7^^'^^ and 07^^^^ + 57^^^"^: 

^5678^±. ^ ^^±. ^ (^^1956 ^ <,^1234)^.± ^ ^^.± ^ (4_-^g) 

Thus the vacuum energy of fermion part is represented by 

Ef = -^{1 + 1 + 3 + 1) = -3fi. (4.19) 

In conclusion, the net vacuum energy is 

^tot = -^l . (4.20) 

Thus the zero-point energy is completely identical with that in the P = case and it is 
independent of the parameter P. 

We can derive the spectrum of the abelian part, and we might expect that the resulting 
spectrum would agree with that in the linearized supergravity in eleven dimensions around 
the time-dependent background consider here (as is shown in case of the eleven-dimensional 
pp-wave background p^]). 



5. Classical Solutions 



In this section we will study the classical solutions of our matrix model. There are the 
extra terms with fij, and so we can expect that graviton solutions would be modified. Also, 
fuzzy sphere solutions should exist due to the presence of the Myers term. 

In order to study classical solutions we need equations of motion. We will set oj = ip 
= for simplicity in the following consideration. By taking a variation of action, the 

^The choice of chirality matrices is not relevant to the projection operator Pi and P2. 
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equations of motion are obtained as 



X' = 


-[[X\X%X'] 


+ 2nPX^ - 2ifi{X^X^ - X^X^) , 


(5.1) 


x^ = 


-[[X^,X'],X'] 


- 2ijlPX^ - V(l - P2)X2 - 2iii{X^X^ - X^X^) 








(5.2) 


x^ = 


-[[X^X%X'] 


- fi'^X^ + 2iixP{X^X^ - X^X^) , 


(5.3) 


x^ = 


-[[X\X%X'] 


- /x^X^ + 2inP{X^X^ - X^X^) , 


(5.4) 


X' = 


-[[X^X%X'] 


- 2ifi{l - P2)V2(J^6J^2 _ j^2j^6) ^ 


(5.5) 


x^ = 


-[[X^X%X'] 


- 2in{l - P2)l/2(j^2j^5 _ j^5j^2) ^ 


(5.6) 


X' = 


-[[X',X%X'] 


-2i,,il-P^y/\X^X^-X^X^), 


(5.7) 


x^ = 


-[[X',X%X'] 


-2ifi{l-PY/\X^X' -X'X^), 


(5.8) 


x^ = 


-[[X^,X%X'] 


- 2ifj,{X^X'^ - X^X'^) + 2ifj,P{X^X'^ - X^X^) . 


(5.9) 



We can evaluate the energy of the classical solution by using the classical potential V in 
our model: 



2V = Tr 



^[X\ X^]^ + V(l - P2)(X2)2 + ;,2(^3^2 ^ ^2(^4^2 



+ ■ siix^'x^x^^ - pxi^x^x^i 

3 

+ (1 - p2)l/2^[2^5^6] ^ (1 _ p2)l/2^[2^7^8l) 



(5.10) 



We will discuss below some types of classical solutions and evaluate the classical ener- 
gies by using the above classical equations of motion and potential. 



5.1 Graviton Solutions 

To begin with, we will consider the graviton solutions. In the case of a matrix model in 
fiat space, the graviton solution is a free particle moving with a constant velocity, while in 
the pp-wave case this solution is modified to be described by an oscillator because of the 
presence of mass terms in the action. In our case the graviton solution is slightly different 
from that in the usual pp-wave cases as we will explain. 

Let us investigate the graviton solution by considering the diagonal matrix. The gravi- 
ton solutions for X^, ■ ■ ■ , X^ are free particles as in flat space. For X^ and X^ the solution 
is a couple of harmonic oscillators as in the usual pp-wave case because of mass terms. The 
remaining parts are X^ and X^. These directions couple each other as follows: 



X^ = 2/xPX2 , X2 ^ _2nPX^ - 4/^2(1 _ p2)x2 . 



(5.11) 



We can easily solve this system of equations and the solution is 



X^ = Px2 sin(2//r) - P^ cos(2//r) + (l - P^) cqt + Cq , 
X^ = X2 cos(2/xr) -I- ^ sin(2/xr) — . 



(5.12) 
(5.13) 
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The X2 and V2 are the initial position and velocity of X'^, respectively. The cq and Cq are 
arbitrary constants related to the initial velocity and position of X^ , respectively. It should 
be remarked that the X^ direction behaves as a harmonic oscillator in spite of the absence 
of its mass term, because of the contribution of the term fij. Notably, the X^ behaves as 
a freely moving particle in the P = case and as a harmonic oscillator in the P = 1 case. 
This result shows a specific feature of our model. 

In addition, we can easily see that there exist a rotating solution: 

X^{t) + iX^T) = {X^{0) + iX\0))e'''^ , [X^{0),X\0)] = , (5.14) 

as in the matrix model on the pp-wave background . The energy of this solution is given 

by 

2V = fi^Ti [{X^{0)f + {X\0)f] . (5.15) 

The above solutions describe the motion of a D-particle in our model. For non-diagonal 
cases, the D-particle can expand to form a higher dimensional object with two dimensions 
(i.e., membrane) because of the Myers terms. The expanded fuzzy membrane solutions 
will be discussed in the next subsection. 

5.2 Fuzzy Sphere and Hyperboloid Solutions 

Here we would like to find fuzzy sphere type solutions by imposing the ansatz: 

X^ = aJ^ , X^ = (3J^ , X^ = aJ'^ , (5.16) 

where the J" (a = l,2,3)'s are the SU{2) generators and satisfy the SU{2) Lie algebra 
[J°-,j'^] = ieabcJ^- By putting the ansatz into the equations of motion, we obtain two 
conditions 

aP'^ + - 2fial3 = , (3 - ^xa^ + 2^m^{1 - P'^)f5 = Q . (5.17) 
These algebraic equations can be easily solved, and we obtain two nontrivial solutions: 

= -2^?{l - p2) + /^(i + ^i + 8(l-P2)), 
/3 = ^^-f\/l + 8(l-P2), (5.18) 

and 

«2 ^ _2^2(i _ p2) ^ /^(i _ ^1 + 8(1-P2)), 

/? = ^^ + f\/l + 8(l-P2), (5.19) 

and a trivial solution a = /? = 0. The first nontrivial solution describes a fuzzy ellipsoidal 

sphere and the second one represents a fuzzy hyperboloidll in the range < P < 1**. In 

"The ellipsoidal sphere and hyperbolic type solutions were constructed in the pp-wave matrix model by 
D. Bak @. 

**If we consider the region P < for the first case, then the shape of the solution becomes a fuzzy 
hyperboloid. On the other hand, if we consider the region P > 1 for the second case, then the solution 
becomes fuzzy ellipsoidal sphere. 
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the second case the value of a becomes purely imaginary and the compact SU{2) group is 
replaced with the noncompact group S0{2, 1). The energies of the fuzzy ellipsoidal sphere 
( |5.18 ) and hyperboloid ( 5.ig| ) are represented by 

2V = (a^/?^ - ^^^a'p^ Tt[{j'f + {j'f] 

+ (a^ + V(i _ p2)^2 _ 4^^2^^ Tr(j2)2 

= -0^13"^ /N, (5.20) 

where we have normalized the trace for the A^-dimensional representation of SU (2) gener- 
ators as 

T¥( JV^) = ^6"^ . (5.21) 

The energy for the fuzzy ellipsoidal sphere ( |5.18| ) is negative while that for the fuzzy 
hyperboloid ( ^.19 ) is positive. 

Other fuzzy solutions are given by 

= c • /3 J\ = c • a J^ 

= c - aJ'^ , otherwise zero , (5.22) 

and 

= c • /3J\ = c • a J2 , 

= c ■ aJ^ , otherwise zero . (5.23) 

In the same way, these two cases lead to the same conditions: 

2c^ {a'^P + 2^2^ - fia"^) = , c^a (2/i/5 - p"^ - a"^) = . (5.24) 

The solutions of the above algebraic equations are 

= -3^2^ /3 = 3fi, (5.25) 

and the trivial solution a = P = 0. Thus, both of the nontrivial solutions (5.22) and ( |5.23| ) 
with (5.25) are fuzzy hyperboloids and have the same energy given by 

2y = ■ + _ 1^«2^^ Tr(ji)2 

= ■ ( + a^- Afia^p + 2a^p'^) /N , 

= 27/c^/iV . (5.26) 

We also can apply the above consideration to the X^, X^ and directions by sup- 
posing the following ansatz: 

X^ = aJ^ , X^ = aJ^ , X^ = PJ^ , otherwise zero . (5.27) 
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By inserting this ansatz into the equations of motion, we can obtain the two constraint 
conditions: 



fia'^P + a'^/3 = 0, + a/?^ + ;u^a + 2/iPa/3 = . (5.28) 
These equations can be solved and the solution is 

a"^ = -^^(l- P^)^ P = _^P_ (5.29) 

This is also fuzzy hyperboloid solution, whose energy is represented by 

2V = l^a'p' + + ^i^Pa'P^ Tr[( Ji)2 + [J^f] + (^a' + ^fiPa'P^ Tv{J^f 

= P^f/N <0. (5.30) 

As the result, the energy of the fuzzy hyperboloid solution has negative energy. 

Now we will discuss the range of the parameter P. For example, if we consider the 

P2 > 1 

case in the solution ( |5.29| ), then the values of a and /3 give the ellipsoidal fuzzy 
sphere solution. If < 1, then the shape of the solution is hyperbolic. This result implies 
that the shape of the solution should be modified as the parameter P changes. It should 
be remarked that the parameter P might run from to 1. If P is out of this range, then 
the bosonic mass term would be tachyonic and it seems that some trouble would appear. 
However, such an issue cannot be caused in the Rosen coordinates since it would be intrinsic 
to the system described in terms of Brinkmann coordinates. 

We will not discuss further classical solutions but it would be expected that other 
interesting solution can be constructed in this model, and so the construction of such 
solutions are an interesting future problem. The stability of the fuzzy hyperboloids, which 
would be non-supersymmetric, would be also interesting^^. 

5.3 BPS Equations 

In the previous subsection, we have presented examples of classical solutions. Now we will 
consider the conditions which lead to supersymmetric classical solutions in terms of BPS 
equations. 

We can derive the BPS equations from the supersymmetric condition: d^ip = 0. The 
conditions preserving all of dynamical supersymmetries are given by 

[X^X^]=i^f^X\ [X^X'] = i^f,X\ [X\X^]=i^f^X\ 

[X\ X'] = -^s ■ ^fsX^ , [X',X^] = -is ■ ^f,X^ , [X^X'] = -^s ■ ^f,X' , 

[X^ , X^] = ic ■ ^a^X^ , [X^ , X^] = ic ■ ^nX^ , [X^ , X^] = ic ■ ^nX^ , 
o o o 

[X'',X^]=ic-^fi{l-a)X\ [X\X^]=ic-^fiX\ [X^,X'']=tc-^fiX\ 
DrX^ = -2nsX^ , DrX^ = 2/isX2 , D^X^ = fiX^ , DrX^ = -^iX'^ , 

DrX^ = DrX^ = DrX'^ = DrX^ = DrX^ = , (5.31) 



^^The quantum stability of fuzzy sphere is discussed in p7[ . 
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where a constant parameter a takes the value in the range < a < 1. The commutators 
which do not appear in the conditions ( 5.31| ) should vanish (for example, = 0). 

We note that the equations ( ^.31 ) are not necessary but sufficient conditions for the BPS 
conditions. Because we are dealing with a non-maximally supersymmetric case, some 
equations are satisfied trivially. 

Now we can find a solution of the above BPS conditions ( 5.31| ) represented by 



X^ 
X^ 
X' 



-/i J-^ ) I (g) ( c • J 



X 



7 



= 101® 
= I(gI(g)I(g) 

{2fi ■ sJ^ 




(5.32) 



The symbol (g means the embedding of the SU{2) generators into the matrix, and the ith. 
generator belongs to the iVj-dimensional representation where = A^i + ^^^^2 + -^3 + -^4- 
For example, (|/iJ^) (g I g) (c • gi (c • |/^</'^) is represented by the following matrix 



c • -/iJ^ (g ( c • -/iJ' 



v 



\ 



>N2 



The above solution ( p.32| ) preserves 4 dynamical supersymmetries while the 16 kinematical 
supersymmetries are broken. Hence this is a l/5(=4/20) BPS object^. In our case this 
configuration has zero energy as in the pp-wave matrix model (i.e., ^ = 0). Therefore, this 
configuration can appear in the classical vacuum without the cost of energy, although this 
is not a classical solution of equation of motion. 



6. Conclusion and Discussion 

We have discussed a matrix model on a homogeneous plane-wave background. The action 
of the matrix model has been proposed. This theory has 20 supersymmetries, and we 
have explicitly constructed the transformation laws of 4 dynamical and 16 kinematical 
supersymmetries in this model. Then, the vacuum energy of the abelian part has been 

■'■■'■In the pp-wave matrix model, such fuzzy sphere solutions preserves 16 dynamical supersymmetries and 
are the l/2(=16/32) BPS objects. 
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calculated, and we have shown that the vacuum energy in the P 7^ is identical with 
that in the P = case. In particular, the effect of a parameter P is to shift the origin of 
the harmonic oscillator and to rescale the continuous spectrum of . We have also 
found classical solutions. The graviton solution is slightly different from that in the usual 
pp-wave cases because of the presence of the additional parameter P. It has been shown 
that the fuzzy hyperboloid and fuzzy sphere solutions are also exist in our model. Notably, 
the fuzzy ellipsoidal sphere solution ( |5.16 ) with ( |5.18| ) and the fuzzy hyperboloid solution 



( 5.27 ) with ( 5.29|) have negative energies. 



It would also be interesting to study a supermembrane theory by taking a large N 



limit (for the pp-wave case, see |9|). By assuming the matrix regularization |1C], the action 
of supermembrane theory can formally be written down as follows: 



S = J dTcPaw{a) - 



-4/i2(l - P'^){X'^f - ^2(^3)2 _ /i2(x4)2 + 2fiPX^DrX'^ - IflPX'^D^X^ 

(M^ - py2) ^ + 2^ . I ^ -p 

I,J,K=1,2,9 I,J.K=3,4:,9 



+ (l-P2)l/2 ^ +(l-p2)l/2 ^ }ejjKX^{X',X'} 



I,J,K=2,5,e I,J,K=2J,8 
^ ^ ^129 _ p^349 ^ _ p2)l/2^256 _^ _ p2)l/2^278 



(6.1) 



where the covariant derivative is Dr* = dr * *} and {A, -B} = {eab/'w)daAdhB is the 
Lie bracket. The indices a,b represent the spacial directions of membrane world- volume. 
An interesting subject to study is the brane charge in our model. The brane charges are 
modified in the nontrivial background [P, |2^. In addition, it would be interesting to study 
an open supermembrane theory on the time-dependent plane- wave background as in 30 1 



because this background has quite special properties and so the classification of Dirichlet 



branes are also nontrivial. All of the above topics will be reported in our next paper |22]. 

One of the most interesting subjects is the stability of a single supermembrane on the 
time-dependent background. In the flat case a single supermembrane is unstable as is well 
known [^]. This fact is deeply based on the cancellation of zero-point energies between 
bosons and fermions. However, it is not clear in general whether the bosonic zero-point 
energy cancels out the fermionic one. Remarkably, our result on the vacuum of the abelian 
part of the matrix model here should support that the zero-point energies cannot cancel 
out. Moreover, the flat directions are left in the time-dependent homogeneous plane-wave 
in contrast with the maximally supersymmetric pp-wave case. This issue is an interesting 
future work. 

Finally, we comment on the relation of our matrix model to the time-dependent back- 
ground as a closing remark. The homogeneous plane-wave background we have considered 
can be mapped to the explicitly time-dependent plane-wave background via the time- 
dependent coordinate transformation. However, once we move to the stationary frame 
in which the time dependence is not manifest, the effect of the time dependence can be 
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expressed by a kind of the vector potential, which leads to the constant magnetic field. 
Thus, we might expect the relationship of our model to the noncommutative geometry or 
Seiberg-Witten map [31|. In addition, the plane- wave geometry considered here is closely 
related to the Lewis- Riesenfeld invariant theory |^2| as noted in This theory leads us 
to the quantization of the system of time-dependent harmonic oscillators. Study of such a 
system would be a useful laboratory to promote the understanding of the time-dependent 
backgrounds and to develop the techniques to treat theories on them. We believe that our 
model should be a clue to shed light on the physics on the time-dependent backgrounds. 
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A. Matrix model on the general homogeneous plane-wave 



Here we shall discuss the matrix model on the general homogeneous plane- wave background 
(anti-Mach space) whose metric is and the four-form flux ( |2.2| ) is equipped. As 

discussed in [^], this background admits extra Killing spinors satisfying 



ex 



(A.2) 



where we have defined as 



(A.3) 



and eo is a constant spinor. 

The matrix theory action on the background is given by 



5=1 IdrTi 



+ ^^|;^Dr^^J - ^|^^J^[X\i>] + l^itl;^{W + + hfi^,,kX'X^ X' 



(A.4) 



where the matrix W and the covariant derivative are defined as, respectively. 



W = — ^ 6 , Dr* = dr * —i[uj, *] ■ 



(A.5) 



The action ( [A.4| ) is invariant under the dynamical supersymmetry transformations 



DrX% - 2f,fi,X'j^ - -[X\X^hij 

-l^xSiw-3cp)-^x\w + ci>H 



(A.6) 
(A.7) 
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provided that 

-36(^,,y + fijfjkl^)]e = 0, (A.8) 
6 = et(^-3'^)-eo. (A.9) 



The Eqs. (^) and ( |Al9|) are nothing but ( |A11 ) and ( |Al2|) respectively. The action 



is also invariant under the kinematical supersymmetries 

5^X^ = V = 0, 5^^^ = ??^, ry(r) = e-t(^+'^No (A.IO) 
where % is a constant spinor. 

B. Analysis of and directions in the P 7^ case 

Here we will examine the and X"^ directions in the P 7^ case where the analysis of 
these directions are quite nontrivial and complicated due to the presence of the nonzero 
fij. The Lagrangian of this part is given by 

L = 1 \^{X^ f + - 4/i2(l - P2)(^2)2 ^ 2flP{X^X^ - X^X^)] . (B.l) 

By using the standard procedure, the Hamiltonian can be derived as follows: 

H = ^{tt^ + fiPXy + ^(tt^ - ^.PX'f + 2f,\l - P^){X^f . (B.2) 

Here we quantize the system by imposing the commutation relation: 

[X\TT^]=i6'^ , 

and will solve the energy eigenvalue problem: 

When we decompose the ^ as \I' = e^*^^^^"''"^^; we can obtain the equation for the x 
described as 

Hx = Ex, 

P ^ i (vri + 2^,PX'f + i(vr2)2 + 2f,\l - P'){X^f . (B.3) 

The expression x = e^^-^^^ (j), where pi is the momentum for the X^-direction, allows us to 
rewrite the Hamiltonian H as 

Thus, we can see that the effect of a parameter P is to shift the origin of harmonic oscillation 
of X^-direction and to rescale the kinetic energy (continuous spectrum) of X^. 

Finally, we would like to note that the final expression of Hamiltonian is closely related 
to the invariant in the context of the Lewis and Riesenfeld theory [^2| , which can quantize 
the system of time-dependent harmonic oscillator. 
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